An improved analytical description of inspiralling and coalescing black-hole binaries 
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We present an analytical formalism, within the Effective-One-Body framework, which predicts 
gravitational-wave signals from inspiralling and coalescing black-hole binaries that agree, within 
numerical errors, with the results of the currently most accurate numerical relativity simulations for 
several different mass ratios. In the equal-mass case, the gravitational wave energy flux predicted by 
our formalism agrees, within numerical errors, with the most accurate numerical-relativity energy 
flux. We think that our formalism opens a realistic possibility of constructing a sufficiently accurate, 
large bank of gravitational wave templates, as needed both for detection and data analysis of (non 
spinning) coalescing binary black holes. 
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The opening of gravitational wave (GW) astronomy 
depends on our theoretical ability at computing, within 
Einstein's theory of General Relativity, a sufficiently ac- 
curate approximation to the GW signal emitted by the 
premier target of ground-based GW detectors: inspi- 
ralling and coalescing binary black holes (BBH's). In- 
deed, to detect these GW signals, and extract physi- 
cal information from them, one must correlate the noisy 
output of the detectors to a very large bank of "GW 
templates" , each template giving an accurate represen- 
tation of the gravitational waveform emitted by a BBH 
with certain physical parameters (notably masses and 
spins). Recent breakthroughs in Numerical Relativity 
(NR) [l], 0, @, 0, H, [f| have given us access to an accu- 
rate knowledge of the waveform emitted during the late 
inspiral and merger of a sparse sample of BBH systems 
(see for a review). However, BBH simulations arc 
time consuming. This precludes the sole use of NR sim- 
ulations for building the needed bank of GW templates, 
densely filling the multidimensional space of BBH phys- 
ical parameters (masses and spins). It is urgent to have 
in hands an analytical method able to give a sufficiently 
accurate representation of the motion of, and the gravita- 
tional radiation from, coalescing binary black holes with 
arbitrary masses and spins. We shall describe here a for- 
malism which hopefully solves this problem, in the case 
of circularized, non-spinning binary black holes with ar- 
bitrary masses mi, m^. 

The analytical formalism presented here is a signif- 
icantly improved version of the general Effective-One- 
Body (EOB) method @, S OS EI- The predictions 
of previous (less accurate) implementations of the EOB 
method have already been compared, with success, to 
various types of results from numerical simulations flS 
03 OS M, OS OS OS OS HI HI- Before explaining 
the improvements that characterize our formalism, let us 
recall that the four essential elements of the EOB ap- 
proach are: (i) a Hamiltonian 7J rea i describing the con- 
servative part of the relative dynamics of the two black 
holes; (ii) a radiation-reaction force T v describing the 



loss of (mechanical) angular momentum, and energy, of 
the binary system; (iii) the definition of the various mul- 
tipolar components of the "inspiral-plus-plunge " ( metric ) 
waveform ^^ pluns °; and (iv) the attachment of a subse- 
quent "ringdown waveform" h"^ down around a certain 
(EOB-determined) "merger time" t m . The latter fourth 
facet of the EOB formalism, namely the assumption of a 
sharp transition, around the BBH merger, between the 
"plunge" and a ringdown behavior, was inspired by the 
classic "plunging test-mass" result of [22j . This assump- 
tion has been well confirmed by the results of NR simula- 
tions @, and we shall not try here to improve on it. We 
shall follow here the usual EOB procedure [9( of matching 
the "insplungc" and "ringdown" waveforms around the 
EOB merger, and matching, time t m , defined as the lo- 
cation of the maximum EOB orbital frequency. [We use 
the specific procedure of pS (3, with 5 quasi-normal- 
modes, and a "comb" of spacing 5 = 2.3Mf , where Mf is 
the mass of the final black hole.] 

On the other hand, we propose here to improve the 
three other basic elements of the EOB formalism in the 
following way: 

(i) The central object entering the relative Hamilto- 
nian, i? rcal = AI[l+2is(H cS -l)} 1 / 2 , where M = m 1 +m 2 , 
(i = mim2/M, v = fJ>/M , z% = 2^(4 — 3^), and 



Heff = ^JpI + A(l/r) (l + p 2 /r 2 + z 3 < /r 2 ) (1) 

is the EOB radial potential A{u) (here u = 1/r; the EOB 
radial coordinate r is rescaled by M, with G = c = 1; and 
p rr is canonically conjugated to the EOB "tortoise-like" 
radial coordinate r*, defined in [3). Current analytical 
calculations, within the post-Newtonian (PN) formalism, 
of the dynamics of BBH's have computed the first four 
terms (3PN approximation) in the (Taylor) expansion of 
the radial potential A(u) in powers of u [lOl. |23||. namely 
A 3PN (u) = 1 - 2u + 2vv? + (94/3 - 4Itt 2 /32) vu A . Here 
we propose to consider the two-parameter class of (ex- 




FIG. 1: Equal-mass case: agreement between NR (black online) and EOB-based (red online) I = m = 2 metric waveforms. 



tensions and) resummations of ^4 3PN (it) denned by 

A(u;a 5 ,a 6 ,v) = P£[A 3PN (u) + va b u 5 + va e u 6 ], (2) 

where P§ denotes a (1,5) Pade approximant. This 
is a generalization of the one-parameter (05) class of 
P4 Pade-resummed A-potentials considered in previous 
EOB works [HI, [T3, E^l - 

(ii) The second (and most novel) improvement that we 
introduce here concerns the radiation reaction force T v . 
We make use of the very recent results of [Hj] concerning 
an "improved resummation" of post-Newtonian multipo- 
lar waveforms. Specifically, we define T v in the following 
way (Jl denoting the EOB orbital frequency): 

^ s -^£I> n ) a i**ffii a - (3) 

1=2 m=l 

Here, we shall take £ max = 8 and define the individ- 
ual multipolar waveforms (where e = 0, 1 labels the 
"even" or "odd" parity) in the following way: (a) the 
leading quadrupolar contribution to T v , i.e. the term 
e = and ^ = m = 2 in (J3J), is computed by using the 
quadrupolar waveform defined in Eq. (QJ below; while 
(b) the subdominant terms (i.e. when, either e = 1, or 
e = and I > 3), are defined by Eq. (1) of [pi, 
together with the other definitions given there (see [25[ 
for details). 

(iii) The third improvement introduced here concerns 
the (dominant) even-parity, quadrupolar (e = and I = 
in = 2) "insplunge" waveform. We take it in the form 

h22 = ^ n 22 xY 2 ^ (|, $) H cS T 22 e*»Mx) fg QC , 

(4) 

where, for notational simplicity, we have suppressed the 
parity label e = 0. Such a multiplicatively decomposed 
form of h 22 was introduced in pjf, [IU . See these refer- 
ences, and [13, H3, [25[ for the definition of the factors 
m EL; f22(x) is a modulus correction, here defined (as 
in [l8[, where f 22 was computed, using [2(| [27|, at the 
3+ 2 p]\f accuracy; Eq. (10) there) as being the following 



Pade-resummed function / P 2 f (a;; v) = P 2 [/ 2 I 2 aylol (x; z^)]. 
The new ingredient of h 22 introduced here is the defi- 
nition of the last factor in Eq. ([?]), namely an additional 
"Next-to-Quasi-Circular" (NQC) correction factor of the 
formdl 

/ 2 N 2 QC ( fll ,a 2 ) = l + a lP 2 r J(rn) 2 +a 2 r/rn 2 . (5) 

A crucial facet of the new analytical formalism presented 
here consists in trying to be as predictive as possible by 
reducing to an absolute minimum the number of "flex- 
ibility parameters" entering our theoretical framework. 
We shall achieve this aim by "analytically" determining 
the two parameters 01,02 entering (via the NQC fac- 
tor Eq. d5J|) the (asymptotic) quadrupolar EOB wave- 
form i?/if 2 OB (where R = R/M) by imposing: (a) that 
the modulus |P/if 2 OB | reaches, at the EOB-detcrmincd 
"merger time" t m , a local maximum, and (b) that the 
value of this maximum EOB modulus is equal to a cer- 
tain (dimensionless) function of v, <p(v). We calibrated 
(independently of the EOB formalism) by extract- 
ing from the best current Numerical Relativity simula- 
tions the maximum value of the modulus of the Numeri- 
cal Relativity quadrupolar metric waveform | P/i^l- Us- 
ing the data reported in Q and [2fj| , and considering the 
"Zcrilli-normalizcd" asymptotic metric waveform ty 22 = 
Rh 22 /V24, we found ip{v) ~ 0.3215^(1 - 0.131(1 - 4i/)). 
Our requirements (a) and (b) impose, for any given A(u) 
potential, two constraints on the two parameters a\,a 2 . 
We can solve these two constraints (by an iteration proce- 
dure) and thereby uniquely determine the values of ai , 02 
corresponding to any given A(u) potential. In particular, 
in the case considered here where A(u; 05, as, v) is defined 
by Eq. this uniquely determines Oi, 02 in function of 

05, ae and v. 

Finally, our analytical formalism contains only two 
analytically undertermined parameters, namely 05 and 

06, which parametrize some flexibility in the Padc- 
resummation of the basic radial potential A(u), con- 
nected to the yet uncalculated higher PN contribu- 
tion [3l|. We have first compared the (05, ag)-dependent 
predictions made by our formalism to the high-accuracy 
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FIG. 2: Phase difference between the analytical and numer- 
ical (metric) waveforms of Fig. Q] 

waveform from an equal-mass BBH (y = 1/4) computed 
by the Caltcch-Cornell group @ (and now made avail- 
able on the web). We found that there is a strong de- 
generacy between 05 and clq in the sense that there is an 
excellent EOB-NR agreement for an extended region in 
the (05, a6)-planc. More precisely, the phase difference 
between the EOB (metric) waveform and the Caltech- 
Cornell one, considered between GW frequencies Mlo-^ = 
0.047 and Mu R = 0.31 (i.e., the last 16 GW cycles be- 
fore merger), stays smaller than 0.02 radians within a 
long and thin banana-like region in the (05, ae)-planc. 
This "good region" approximately extends between the 
points (05, a 6 ) = (0, —20) and (as, qq) = (—36, +520). As 
an example (which actually lies on the boundary of the 
"good region" ) , we shall consider here the specific values 
as = 0,a 6 = —20 (to which correspond, when v = 1/4, 
ai = -0.036347, a 2 = 1.2468). We henceforth use M as 
time unit. 

Figure[T]compares (the real part of) our analytical met- 
ric quadrupolar waveform , I'f 2 OB /i/ to the corresponding 
(Caltech-Cornell) NR metric waveform iff™ l v (obtained 
by a double time-integration, a la [20l |. from the origi- 
nal NR curvature waveform ■0 22 ). [We used the "two- 
frequency pinching technique" of [l9| with u)\ = 0.047 
and L02 = 0.31.] The agreement between the analyti- 
cal prediction and the NR result is striking, even around 
the merger (see the close-up on the right). The phas- 
ing agreement is excellent over the full time span of the 
simulation (which covers 32 cycles of inspiral and about 
6 cycles of ringdown), while the modulus agreement is 
excellent over the full span, apart from two cycles after 
merger where one can notice a difference. A more quan- 
titative assessment of the phase agreement is given in 
Fig. [2 which plots the (wi-W2-pinched) phase difference 
Acj) = </>m?tric — ^metric- A</> remains remarkably small 
(~ ±0.02 radians) during the entire inspiral and plunge 
(1V2 = 0.31 being quite near the merger, see inset). By 
comparison, the root-sum of the various numerical errors 
on the phase (numerical truncation, outer boundary, ex- 
trapolation to infinity) is about 0.023 radians during the 



FIG. 3: Equal mass case: metric-amplitudes comparison. 
The maximum of the orbital frequency f2 defines the EOB 
merger. 

inspiral @. At the merger, and during the ringdown, 
A</> takes somewhat larger values (~ ±0.1 radians), but 
it oscillates around zero, so that, on average, it stays 
very well in phase with the NR waveform (as is clear 
on Fig. [!}. By comparison, we note that [6j mentions 
that the phase error linked to the extrapolation to infin- 
ity doubles during ringdown. We also found that the NR 
signal after merger is contaminated by unphysical oscilla- 
tions. We then note that the total "two-sigma" NR error 
level estimated in [f| rises to 0.05 radians during ring- 
down, which is comparable to the EOB-NR phase dis- 
agreement. Figure [3] compares the analytical and numer- 
ical metric moduli, | v E , 22|/i / - Again our (Pade-resummed, 
NQC-corrected) analytical waveform yields a remarkably 
accurate description of the inspiral NR waveform. During 
the early inspiral the fractional agreement between the 
moduli is at the 3 x 10 -3 level; as late as time t = 3900, 
which corresponds to 1.5 GW cycles before merger, the 
agreement is better than 1 x 10~ 3 . The discrepancy be- 
tween the two moduli starts being visible only just before 
and just after merger (where it remains at the 2.5 x 10~ 2 
level). This very nice agreement should be compared 
with the previously considered EOB waveforms (which 
had a more primitive NQC factor, with a 2 = [ijj l20j ) 
which led to large moduli disagreements (~ 20%, see 
Fig. 9 in [2(|) at merger. By contrast, the present moduli 
disagreement is comparable to the estimated NR modu- 
lus fractional error (whose two-sigma level is 2.2 x 10~ 2 
after merger 0]). 

We also explored another aspect of the physical sound- 
ness of our analytical formalism: the triple comparison 
between (i) the NR GW energy flux at infinity (which 
was computed in [2l|): (ii) the corresponding analyti- 
cally predicted GW energy flux at infinity (computed by 
summing \hi m \ 2 over t,m ); and (hi) (minus) the me- 
chanical energy loss of the system, as predicted by the 
general EOB formalism, i.e. the "work" done by the ra- 
diation reaction ^mechanical = ^-7>- This comparison is 
shown in Fig. [H which should be compared to Fig. 9 
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FIG. 4: Triple comparison between NR and EOB GW energy 
fluxes and the EOB mechanical energy loss. 
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FIG. 5: Unequal mass case: Comparison between metric 
waveforms for the 2:1 mass ratio. 



of [2l[. We kept the same vertical scale as [2l[ which 
compared the NR flux to older versions of (resummed 
and non-resummed) analytical fluxes and needed such a 
±10% vertical scale to accomodate all the models they 
considered. [The horizontal axis is the frequency w of 
the differentiated metric waveform ft-22-] By contrast, 
we see again the striking closeness (at the ~ 2 x 10 -3 
level) between the analytical and NR GW fluxes. As 
both fluxes include higher multipoles than the (2, 2) one, 
this closeness is a further test of the agreement between 
our analytical formalism and NR results. [We think that 
the ~ 2a difference between the (coinciding) analytical 
curves and the NR one on the left of the Figure is due 
to uncertainties in the flux computation of [211 ] . possibly 
related to their method of computing h.} Note that the 
rather close agreement between the analytical energy flux 
and the mechanical energy loss during late inspiral is not 
required by physics (because of the well-known "Schott 



term" 28J), but is rather an indication that he m can be 
well approximated by —imClhi m (used in Eq. ([3])). 

Finally, as the power of our formalism resides in the 
ease with which it can accomodate continuous varia- 
tions in the basic physical parameters of the consid- 
ered BBH, we shall discuss an unequal-mass system 
(y < 1/4). The highest-accuracy data that we had in 
hands is the Jena-group simulation of a 2:1 mass ratio 
BBH (y = 2/9 = 0.22222). When a 5 = 0, a 6 = -20, and 
v = 2/9, one finds a x = -0.017017, a 2 = 1.1906. Using 
the data reported in [2(| (and the function ip(y) quoted 
above), we compare in Fig. [5] (the real part of) our ana- 
lytical metric quadrupolar waveform ^ , f 2 OB / i/ to the cor- 
responding (Jena, 2:1 mass ratio) NR metric waveform 
^ r 22 ? 'A / - [We use, as in [2(j, the pinching frequencies 
uj\ = 0.1005, W2 = 0.4542 .] Again we have an excel- 
lent analytical- numerical agreement, both in phase and 
in modulus. The small differences between the two are 
within the numerical errors (see (20j ). 

Conclusions. We have described a specific analyti- 
cal formalism (within the EOB framework), which con- 
tains as arbitrariness only the resummation-flcxibility pa- 
rameters of the crucial EOB A(u) potential. We have 
shown that for certain values [321 ] of these parameters 
(as, ae): (i) the waveform predicted by our analytical for- 
malism agrees, essentially within numerical errors, with 
the currently most accurate numerical relativity simu- 
lations; this agreement holds for several different mass 
ratios (1:1, 2:1 and 4:1 - not shown here); and (ii) the 
gravitational wave energy flux predicted by our formal- 
ism agrees, within numerical errors, with the most accu- 
rate numerical-relativity energy flux. We think that our 
formalism (possibly after some further minor improve- 
ments) opens a realistic possibility of constructing (with 
minimal computational resources) a very accurate, large 
bank of gravitational wave templates, thereby helping 
in both detecting and analyzing the signals emitted by 
inspiralling and coalescing binary black holes. [Though 
we have in mind essentially ground-based detectors, we 
think our method could also apply to space-based ones.] 
Finally, from a theoretical point of view, we think that 
our method can be extended to the description of (nearly 
circularized) spinning black hole systems (see [Hj]). 

After the submission of this work, a paper [29| compar- 
ing Caltech-Cornell numerical data to a different version 
of the EOB formalism appeared on the archives. The 
EOB formalism of [29| does not use our novel (predic- 
tive) radiation reaction Eq. ([3]) (but rather the u po ie- 
tuncd one advocated in [l8|), nor our ag-improved A 
potential, Eq. @. Moreover, by contrast to the ap- 
proach advocated here to reduce to an absolute minimum 
the number of adjusted parameters, namely two, (as, 
as), Ref. [2!| tunes six parameters: (05(1/4), v po i e (l/4), 
a 3 ^(l/4), aJ»(l/4), 0^(1/4), and ^ 2 atch (l/4)). 
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Note that one could also similarly improve the subleading 
higher- multipolar-order contributions to T^. We leave 
the exploration of such refinements to future work. 
Indeed, as and formally parametrize the 4PN and 5PN 
contributions to the Taylor expansion of A(u). 
We leave a detailed discussion of the "degenerate set" of 
"good" (a5,ag) values to a future publication [25( [ . 



